Depinning and critical current characteristics of topologically defected vortex lattices 
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We discuss the role of dislocation assemblies such as grain boundaries in the dynamic response of 
a driven vortex lattice. We simulate the depinning of a field-cooled vortex polycrystal and observe 
a general enhancement of the critical current as well as a distinct crossover in the characterisitic 
of this quantity as a function of pinning density. The results agree with analytical predictions for 
grain boundary depinning. The dynamics of grain boundaries thus proves an essential mechanism 
underlying the flow response of defected vortex lattices and the corresponding transport properties 
of the superconducting material. We emphasize the connection between the topological rearrange- 
ments of the lattice and its threshold dynamics. Our theory encompasses a variety of experimental 
observations in vortex matter as well as in colloidal crystals. 



The possibility of tuning mechanical properties of vor- 
tex lattices in type II superconductors has made of vor- 
tex physics a versatile framework to study several central 
problems of condensed matter flow [l| . One of the most 
intriguing experimental features of vortex dynamics in 
disordered samples is the emergence of memory effects. 
Hysteresis of the I — V curve is often encountered when 
driving defected vortex lattices above and below their de- 
pinning threshold in both low Tc superconductors, such 
as NbSe2 0, H, Q and high-T(^nisotropic superconduc- 
tors, such as B2Sr2CaCu208 @. Numerical evidence of 
hysteretic behavior has also been collected, by simulating 
driven vortex dynamics in the presence of quenched disor- 
der. The system was initially "prepared" in a disordered 
state by letting it relax from a high-temperature liquid- 
like phase, according to what is know as a field-cooling 
protocol in experiments @. Afterwards currents were 
ramped up and down and two distinct branches could be 
found in the I — V characteristic. The emergence of hys- 
teresis indeed represents a shattering evidence of plastic 
flow in driven lattices. The voltage recorded in exper- 
iments is a direct measure of collective velocity of the 
vortex ensemble, while the applied currents act on vor- 
tices as an external driving force through a Lorentz-like 
coupling [ij . An I — V measurement is thus an exper- 
imental visualization of the force vs. velocity relation 
in a driven over-damped medium. The elastic theory of 
transport in such systems, the well- know elastic depin- 
ning theory, does not account for history dependence in 
force- velocity relations [7]- Such beyond-elastic features 
suggest the appearance of plastic phenomena and their 
involvement in memory effects. The critical current of 
the defected vortex array is thus a plastic threshold to 
vortex motion, which is found to depend on the history 
of the sample. Several experimental studies of the crit- 
ical current Jc of vortex lattices have been carried out 
through the past decades. They include measurements 
of Jc as a function of both the applied magnetic field 
[3, Q and the temperature ioj. In all systems, as a 
field-cooled vortex assembly is driven by a DC or pulsed 
signal, it behaves as if it was driven out of a metastable 
state through an annealing process, so that the critical 



current of the annealed sample proves lower than that of 
the initial field-cooled configuration. Such a process is 
governed by a current-dependent annealing time, which 
diverges as the critical current is approached from above 
,2]. In fact it was shown that in the presence of surface 
barriers at the edge of a sample, the field-cooled and the 
annealed phases can even coexist, as a result of the bal- 
ance between bulk currents, which drive the annealing 
process, and edge currents, which produce contamina- 
tion effects [l^ [ll| . It is now widely accepted that such 
a non-annealed field-cooled state should correspond to a 
peculiar topological rearrangement of the vortex array, 
namely a disordered phase. This is in agreement with 
the experimental observation of glassy features such as 
power-law distributed relaxation times [?]. High-drive 
annealing is thus a process by which the vortex assem- 
bly recovers its topological order. The correspondence 
between topological disorder and higher critical currents 
is ubiquitous in such systems, another example being the 
phenomenology of the so called peak effect. In that case 
a critical current jump accompanies a disordering transi- 
tion at high temperature or field. In field-cooling exper- 
iments, instead, one focuses on the metastable low-field 
and low-temperature state originated after a rapid tem- 
perature quench. 

As for the nature of the field-cooled phase, Delaunay 
triangulation patterns suggest that the topological order 
of the vortex lattice is broken by edge dislocations, which 
tend to arrange themselves into linear arrays such as 
grain boundaries (GBs). The vortex array is thus frozen 
in a polycrystalline state. Vortex polycrystals have been 
observed, after field- cooliw, in various superconducting 
materials such as NbMo fTl, [11, NbSe2 [ll, [H, [H, [13, 
BSSCO [3 and YBCO jH. 

In this Letter, we address the problem of plastic depin- 
ning of such systems. We simulate driven vortex poly- 
crystals and demonstrate that for weak pinning forces, 
the threshold behavior observed in field-cooled samples 
agrees with our analytical predictions based on mech- 
anisms of grain boundary depinning. We are able of 
establishing a tight connection between the topology of 
the vortex polycrystal and its electrodynamical response. 
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This correspondence was first suggested in early numer- 
ical studies [20I . [2l| . Our aim is to determine how the 
current response of these systems is affected by tuning 
crucial parameters such as the defect density Np and the 
magnetic field (proportional to the number of vortices 
Ny), and to emphasize how the topology of the vortex 
ensemble refiects these changes. Transport properties in 
our simulations are quantified by looking at the critical 
current Jc of the vortex array and the collective velocity 
of the system in the steady state v (more precisely, its 
component along the vortex drift direction). 

Simulations are performed in a system of linear size 
L = 36A [12] with periodic boundary conditions. All 
lengths are expressed in units of A. In particular, we 
choose a value of ^ = 0.2A for the coherence length 
(Ginzburg-Landau parameter found e.g. in 

certain low- Tc superconducting alloys). Vortex dynam- 
ics is simulated by integrating Langevin equations of mo- 
tions in the presence of quenched disorder, following the 
same procedure described in a previous publication 0. 
Field-cooled (FC) configurations are obtained by letting 
a random vortex array relax in the impure environment 
and in the absence of external forces. The systems relaxes 
into a vortex polycrystal, in agreement with experimental 
evidence. 

We expect the number of defects Np to affect the re- 
laxation process. As a matter of fact, defects are respon- 
sible for the polycrystalline order observed after relax- 
ation. Vortex mutual interactions try to restore lattice 
order, however defects hinder this process. While varying 
the number of pinning centers Np, we observed that all 
relaxed systems exhibit grain structure. However typical 
grain sizes are found to decrease while increasing Np only 
for low defect densities. Further increases of Np do not 
produce any appreciable drop in grain sizes. At very high 
Np, grain structure becomes extremely complex, however 
it seems that grains have reached a minimum limiting 
size. Figures [T] (a) and (b) display typical vortex array 
topologies in both regimes. 

We found that this behavior shows a natural corre- 
spondence with critical currents. Starting from the above 
relaxed configurations, we determine critical currents, as 
reported in Figure [T] (c). As expected, we see critical 
currents increase with increasing defect densities. At the 
same time, we observe that where we previously spotted 
a qualitative change in the grain size decrease, a crossover 
from a linear to a slower increase in critical currents ap- 
pears. Similar behavior is obtained for different vortex 
densities. Indeed, Figure O shows the same critical cur- 
rent traits for simulations with different numbers of vor- 
tices in the system N^ = 1020, 2016, and 3120. 

In order to understand this correspondence, we have 
to focus on how grain structure might affect depinning 
of the vortex assembly. Compared to a perfect lattice, 
a polycrystal is a system with a much larger number of 
degrees of freedom, coming from the countless ways dis- 
locations and grain boundaries can cooperatively move 
and rearrange. As a consequence, a vortex polycrystal 
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FIG. 1: (a) Grain structure of the relaxed state for low pin- 
ning densities (A'^p — 1032) and (b) for high pinning densities 
{Np = 8256). Dislocations are identified as pairs of 7— and 
5— coordinated vortices (blue and red respectively - color on- 
line) (c) Critical current of the field-cooled vortex array as a 
function of the number of quenched impurities, as obtained 
for a vortex polycrystal (triangles) and a vortex single crystal 
(diamonds), for N^ = 3120 vortices. Numerical results are 
compared to theoretical predictions (dotted and solid lines) 
for grain boundary pinning (see text). The straight line join- 
ing vortex-crystal points (dash-dotted) is drawn as a guide to 
the eye [ll]. 



will find it easier to adjust to disorder than a perfect lat- 
tice or even than a less defected polycrystal. A better 
adjustment to disorder means in general a higher depin- 
ning threshold and eventually a higher critical current. 
Then it is no coincidence that critical currents increase 
rapidly (linearly) in the region where grain sizes decrease 
and slow down where no further grain size decrease can 
be discerned. 

At this point, one might wonder what the behavior of 
a perfect crystal would be under the same conditions. A 
perfect crystal relaxes in the presence of weak disorder 
into the well known Bragg glass phase, which retains tri- 
angular lattice ordering over accessible length scales 0| . 
Such an experimental protocol is known as zero field- 
cooling (ZFC). We simulated ZFC lattices and measured 
critical currents as we did for their FC polycrystalline 
counterparts. Results for the case of 7V„ = 3120 are also 
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FIG. 2: Rescaled critical current as a function of the number 
of defects for three different vortex densities. The cases of 
A''„ = 1020, 2016, and 3120 vortices are depicted. Numeri- 
cal data collapse where an individual GB pinning description 
holds, while they deviate as soon as collective GB pinning 
mechanisms take over. The dotted line is drawn as a guide to 
the eye. 



FIG. 3: Velocity fields at incipient depinning for a FC poly- 
crystal (left) and a ZFC perfect lattice (right). Plastic depin- 
ning and dislocation dynamics are signaled by the emergence 
of heterogeneity and convective vortex rearrangements. In a 
perfect vortex lattice, instead, elastic depinning is accompa- 
nied by widespread avalanches. 



shown in Figure [T] (c) for comparison. Two main ob- 
servations can be made: i) as discussed above and con- 
firmed by experiments [1, H, 0, Q , critical currents are 
higher for a vortex polycrystal; ii) no crossover behavior 
is observed for a vortex crystal as the number of pinning 
points Np is increased. The latter observation represents 
a crucial result of our work. The dynamic response of 
a vortex polycrystal proves radically different from that 
of a dislocation-free lattice. The presence of dislocation 
assemblies such as grain boundaries seems to affect the 
depinning of the vortex array in dramatic fashion. 

We found that close to depinning the dynamics are 
highly heterogeneous (see Figure [3]) . Several swirls ap- 
pear and disappear intermittently in the vortex veloc- 
ity field and produce significant transversal excursions. 
By superimposing vortex velocity fields and the corre- 
sponding topological arrangements (Figure 2]) and moni- 




FIG. 4: Velocity fields at incipient depinning for a FC poly- 
crystal. The grain structure is highlighted. Vortex motion is 
nucleated exclusively around depinning grain boundaries. 



toring the time evolution of the system, we can conclude 
that vortex motion is activated exclusively around de- 
pinned grain boundaries and complicated velocity pat- 
terns emerge in response to GB and dislocation motion. 
Such regions show higher mobility and large transver- 
sal deviations. On the contrary, dislocation free regions 
remain dynamically frozen, unless forces well above the 
depinning threshold are applied. Heterogeneity and co- 
existence of pinned and moving regions are indeed ob- 
served experimentally in vortex matter (25j as well as in 
colloidal polycrystals [11] . 

Such observations corroborate the idea that grain 
boundaries and possibly dislocation assemblies in general 
are crucial in the depinning transition of vortex poly- 
crystals. As a consequence, we propose that the thresh- 
old mechanisms of a driven vortex polycrystal should be 
governed by grain boundary depinning. In two recent 
publications of ours we proposed a theory of GB depin- 
ning [6„,27J. For high defect concentrations, the vortex 
lattice is collectively deformed by the superposition of 
strain fields produced by defects. The straining of the 
vortex lattice results in a collective stress field acting on 
vortex lattice dislocations and grain boundaries as an ef- 
fective pinning field. We refer to this regime as collective 
GB pinning and recall that the critical stress ctc for GB 
depinning under the conditions of our simulations obeys 
the following relation ^6|| 
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(1) 



where D is the average dislocation spacing within the 
GB and b the modulus of the dislocation Burgers vector. 
Ra is the collective pinning length of a 2D vortex lattice 

1/2 1/2 i L 

and its expression is given by Ra — cg()/{fonp nj ) [20|, 
where /o is the typical pinning strength acting on the lat- 
tice, rip is the defect density, strictly proportional to A'p, 
flu is the vortex density, while cge is the shear modulus 
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of the vortex lattice. K is instead the shear modulus of 
the vortex polycrystal. Under such conditions of high 
pinning density, we already pointed out that the vortex 
lattice is highly defected and one has to assume that in 
principle K ^ cge- For low defect concentrations, in- 
stead, pinning is not mediated by the lattice and grain 
boundary dislocations are individually pinned by impu- 
rities. The individual GB depinning stress then reads [2 7| 

K' is the shear modulus of the vortex polycrystal for 
low pinning concentrations. Under such conditions, the 
vortex lattices is split in very large grains and one can 
assume that, to a first approximation, K' « cqq. Re- 
placing the correct expression for Ra and focusing on the 
dependence on Np, Ny and /o we can conclude that for 
the critical current of a vortex polycrystal, the following 
relations hold 

/ /o {^p/^v) individual GB pinning 

\ /o {Np/Nyy/^K{N„)/Ny collective GB pinning 

respectively for low and high defect concentrations. In 
Equations ([3]) we have exploited the proportionality of 
Ny, cge and the magnetic induction B. It is also im- 
plicit in ([3]) that the critical current is proportional 
to the 2D critical force bac- We can now compare our 
theoretical predictions of grain boundary pinning to our 
numerical results for FC driven vortex polycrystals. Fig- 
ure [T] shows how, for low defect concentrations, critical 
currents grow linearly with Np, while for higher pinning 
densities Jc crosses over to a square-root growth. The 
agreement with the theory is remarkable. At this point, 
we should remark that in our simulations grain bound- 
aries develop spontaneously during the relaxation pro- 
cess. The fact that the global behavior matches the one 
analytically predicted for a GB-pinning dominated sys- 
tem allows us to conclude that GB pinning is the relevant 
mechanism that drives vortex polycrystal depinning and 
produces critical current anomalies, such those observed 
in experiments. The first equation in ([3]) also allows a 
prediction of the dependence of Jc on the magnetic field 
H oc Ny , for the case of low pinning density. Our theory 
suggests Jc oc 1/Ny. Indeed, Figure [5] shows that our 
numerical results appear in good agreement with that 
prediction. Experimental curves of Jc vs. H for FC sam- 
ples qualitatively confirm that behavior, for relatively low 
fields, far enough from the peak-effect region [5, 's']. How- 
ever to our knowledge, no quantitative data are available 
at present. This would represent an interesting experi- 
mental test of our prediction, together with a study of 
the Jc vs. H relation in the collective pinning regime. 



This would unveil the physics hidden behind the K{Ny) 
function in Equation ([3]) and ultimately allow a deeper 
understanding of how a complex GB network statistically 
affects mechanical properties of a vortex lattice. 

In the light of the above observations, we can infer 
that plastic flow in vortex polycrystals is dominated by 
the gliding motion of dislocation assemblies such as grain 
boundaries, in analogy with nanocystalline materials. A 
typical fingerprint that corroborates our view is the ex- 
perimental observation of broad-band noise spectra in 
DC-driven vortex arrays^]. Such a feature, also known 
as 1// noise, is a natural consequence of the collective 
dynamics of dislocation assemblies [29j, as opposed to the 
washboard- frequency peaks, which commonly encompass 
the periodicity of a dislocation-free lattice. Indeed, by 
means of Delaunay triangulations we observe that in our 
simulations, right above depinning, dislocations are cre- 
ated and annihilated at the same rate. At higher drives, 
in agreement with experiments, healing takes place in 
the form of a recrystallization process. After healing, 
by lowering the applied current the hysteresis cycle is 
(3) recovered [6|. Current annealing is more effective for 
higher vortex densities, i.e. stiff er vortex arrays, hence 
the steady flow of a dense polycrystalline vortex array 
along the driving direction eventually resembles the flow 
response of a perfect vortex crystal. On the contrary, less 
dense or softer arrays remain topologically defected for 
all the timespan of the numerical simulations. 

In conclusion we performed a numerical study of plas- 
tic depinning in field-cooled vortex lattices. We observed 
the emergence of dislocation assemblies such as grain 
boundaries and investigated their role in vortex dynam- 
ics. Critical currents are found to follow the laws that 
we predicted analytically for grain boundary depinning 
processes. Grain boundary depinning thus proves the 
relevant mechanism accounting for the response of such 
systems. Our theory of GB depinning is able to explain 
several experimental observations like the spatial hetero- 
geneity of the system at the threshold, the high-drive 
annealing of the FC state, the emergence of memory ef- 
fects and the appearance of broad-band voltage noise. 
We established a tight connection between the topology 
of the vortex polycrystal and its response to externally 
applied currents and explored the dependance of critical 
currents and velocities on the density of impurities. Re- 
cent experiments on colloidal polycrystals have partially 
studied this dependance and provided results that agree 
with the predictions of our model [26j . We would like to 
propose a more systematic analysis in the same direction. 
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